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Abstract 

We consider the reactive Boussinesq equations in a slanted cylinder, with zero stress boundary 
conditions and arbitrary Rayleigh number. We show that the equations have non-planar traveling 
front solutions that propagate at a constant speed. We also establish uniform upper bounds on 
the burning rate and the flow velocity for general front-like initial data for the Cauchy problem. 

1 Introduction 

The existence of traveling fronts for reaction-diffusion equations and their stability has been exten- 
sively studied since the pioneering work of Kolmogorov, Petrovskii and Piskunov |26| and Fisher 
|17j . A large number of results have been obtained during the last decade on the generaliza- 
tion of the notion of a traveling front to reaction-diffusion-advection equations in a prescribed 
flow. These include non-planar traveling fronts in shear flows I1D| I12j . and pulsating travel- 
ing fronts in periodic flows [ffl 1391 HTJ] . as well as results for monotonic systems in a unidirectional 
flow jH3 IM1 E3 ■ One of the main qualitative effects of a flow is the speed-up of front propagation 
due to front stretching. Various bounds have been obtained for the speed of propagation of fronts 
in prescribed flows [IJ [21 El 13 El HH 123 HH1 [23 ED , including variational principles for the front 
speed [HI El El EU [22] • The homogenization limit in a periodic flow has also been studied [2"7| . 
Extensive recent overviews can be found in [51 I3U1 l4"T] . 

However, those results have been obtained under the assumption that the flow is imposed from 
outside, and that it is not affected by the evolution of the solution of the reaction-diffusion-advection 
equation, that is, by the temperature or concentration of the reactant. This is known as the constant 
density approximation in the combustion literature. A first step in the coupling of the temperature 
and fluid flow evolution is via the Boussinesq approximation: the density mismatch is so small that 
the density difference is accounted by a buoyancy force in the equation for an incompressible flow. 
Recently a number of works considered systems of a reaction-diffusion-advection equation coupled to 
a flow equation of the Boussinesq type. Global existence and regularity of solutions in two dimensions 
was studied in [2H]- It has been shown that non-planar convective traveling fronts may not exist in 
a vertical cylinder if the Rayleigh number is too small while for large Rayleigh numbers the planar 
fronts become unstable |141 1321 l3*3*j . Moreover, there exists a bifurcation at a critical value p c > 
- non-trivial convective fronts may exist for the Rayleigh numbers close to p c |32| I33j . Numerical 
computations |34| show that non-planar convective fronts exist and are stable for a large range of 
Rayleigh numbers p > p c . The fingering instability in this regime was investigated in |15j . 
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One of the difficulties in the analysis of the Boussinesq problem at large Rayleigh numbers in a 
vertical cylinder is the presence of unstable planar fronts that make uniform lower bounds on the 
front speed quite difficult. However, it has been observed in 0j that such planar fronts cannot exist 
in a horizontal cylinder. One of the main results of [I] is that non-planar fronts in a horizontal 
cylinder exist for small Rayleigh numbers. A purpose of the present paper is to extend this result 
to all positive Rayleigh numbers; we use an approach that is different from [1] and is based on the 
a priori bounds developed in |14j . 

The reactive Boussinesq equations for the temperature T and flow u have the dimensional form 

v 2 

It + u • VT = kAT + —f(T) (1.1) 
u t + u • Vu — uAu + Vp = gTe z 

V • u = 0. 

Here e 2 is the unit vector in the vertical direction, g is the strength of gravity, the speed vq is 
proportional to the traveling front speed in the absence of gravity, k is the thermal diffusivity and 
v is the fluid viscosity. The temperature is normalized so that < T < 1. The nonlinearity /(T) is 
assumed to be a Lipschitz function of the ignition type 

f(T) = for < T < 9 with 9 > 0, f(T) > for T G (0 , 1) and /(l) = 0. (1.2) 

We consider the equations Q1-1JI in a slanted two-dimensional cylinder ax < z < ax + H 

with a finite slope a < oo. It is convenient to rotate the cylinder in order to make it horizontal to 
simplify the notation. Then Ql.l|) becomes 

v 2 

T t + u- VT = kAT + —f(T) (1.3) 

K 

ut + u • Vu — vAu + Vp = gTe 

V • u = 0, 

where u is the flow velocity measured relative to the new coordinate system. The gravity on the 
right points in a direction e that is non-parallel to the the original cylinder was assumed 

to be non-vertical (a < oo). The new rotated problem is posed in a cylinder D = M x x [0,L] 2 , 
L = H/y/l + a 2 . The boundary conditions for the temperature T are set to be front-like: 

dT 

T — > 1 as x — > — oo, T — > as x — > +oo, -7— = at z = 0, L. (1-4) 

oz 

The flow u = (v, w) satisfies the no stress boundary conditions: 

u,u->0as3;-> ±00 and w, uj = at z = 0, L. (1-5) 
Here u = w x — v z is the flow vorticity so that 

Av = —lj z , Aw = UJ X . 

In order to pass to the non-dimensional variables we introduce the laminar front width 5 = k/vq and 
reaction time t c = k/vq and rescale the space and time variables: x new = and t new = t id/t c . 

We also rescale the flow u new = u ^/?;o. Then the Boussinesq equations become 

T t + u • VT = AT + f(T) (1.6) 
Ut + u • Vu — <tAu + Vp = pTe 

V • u = 0, 
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where a = v/k is the Prandtl number and p = gd 3 / n 2 is the Rayleigh number. The problem is now 
posed in the strip D = M x x [0, X] z , A = L/5, with the boundary conditions that come from Q1.4JI 
and (fT3|) . 

The traveling front solutions of 1)1.6)) are solutions of the form T{x — ct, z), u(x — ct, z) with the 
speed c to be determined. They satisfy 



with the boundary conditions 



— cT x + u • VT = AT + f(T) 
—cu x + u • Vu — erAu + Vp = pTe 

V • u = 0, 



(1.7) 



T — > T_ as x — > —oo, T — » as x — > +oo, — — = at z = 0, A (1-8) 

oz 

and 

«),w = 0atz = 0,A. (1.9) 

Here T_ is a constant that is not a priori prescribed. We recall that, as has been observed in [I], if 
the direction of gravity e is not parallel to the x-axis, any traveling front solution of 1)1.7)) must be 
non-planar, that is, it must depend on both variables x and z. This is the main difference between 
the cases of a vertical and slanted cylinder: planar fronts exist in the former case but not in the 
latter. Our main result is the following theorem. 

Theorem 1.1 Let the nonlinearity f(T) be of the ignition type Then a traveling front solution 

(c, T, u) of jl.T\ ) exists such that it is non-planar: T z ^ 0, the flow u ^ and the reaction rate 
f(T) ^ 0. Moreover, the solution satisfies the following properties: c > 0, T G C 2 ' a (D), VT £ 
L 2 (D), u e H l (D) n C 2 ' a (D). If we assume in addition that 

f(T)<(T-6 ) 2 + /X 2 , (1.10) 

then the left limit T_ = 1. 

The assumption (|1.1U|) is of technical nature. It does not involve the Rayleigh number p, it is rather a 
restriction on the channel width A. We do not address the question of the uniqueness of the traveling 
front speed or profile in this paper - this problem requires an additional study. Our results can be 
generalized to the no-slip boundary conditions u = on dD at the expense of a more technical proof 
- we leave this problem for a future publication. 

The general idea of the proof is as follows. We first consider the problem (|1.7|) on a finite domain 
D a = [—a, a] x x [0, A]. Solutions (T^, uJJ) of the restricted problem exist for all c£l. We normalize 
them by the requirement that 

maxT^(x,z) = e . (1.11) 

x>0 

This imposes a restriction on the speed c. In order to show that there exists a speed c a so that 1)1.11(1 
holds we first obtain some a priori bounds on c, T and u under the condition ()l.ll(l . Then we use 
the Leray-Schauder topological degree theory and the above a priori bounds to show that c a exists. 
The a priori bounds allow us to pass to the limit a — > oo. Finally we show that the right limit of T 
as x — > +oo is equal to zero, and that the left limit is equal to one under the additional assumption 
on f{T) in Theorem ll.il This general strategy is similar to that in the proof of existence of traveling 
fronts in a prescribed decoupled flow, as in, for example, [9*lll2j. The main difficulty and novelty are 
in the a priori bounds for the solution of the coupled problem in a bounded domain. 
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Our second result shows that the solution of the Cauchy problem for (|1.6|) propagate with a finite 
speed and that this speed is close to the speed of the laminar front cq when the Rayleigh number is 
small. Recall that there exists a unique speed Co so that a traveling front solution of 

-c^ x = <5> xx + $([/), $(-oo) = 1, $(+oo) = 

exists. 

In order to make this precise we define the bulk burning rate V(t), the Nusselt number N(t) and 
the average horizontal flow U(t) by 

V(t) = 7 T I V ( s ) ds > V(t) = [ f(T)^, (1.12) 



tj J w ' w J A 
N(t) =\f Q N(s)ds, N(t) = J | VT| 2 ^, (1.13) 

i r* 

U(t) = - / ||«(s)||oorfs. (1-14) 
* io 

The following theorem provides uniform bounds on these bulk quantities. It also shows that the 
coupled problem ()1.6|) is in a sense a " regular perturbation" of the single reaction-diffusion equation 
with p = 0. 

Theorem 1.2 Assume that there exists R so that Tq(x,z) = for x > R and Tq(x,z) = 1 for 
x < —R and that the initial vorticity ujq G L 2 (D). There exists a constant C > so that under the 
above assumptions on the initial data Tq, uq we have the following bounds 

c - C[p + p 2 } + o(l) < V{t) < c + C[p + p 2 ] + o(l) (1.15) 

2 



<"/>+ v/| + CV 



+0(1) 



i\T(t)< 

C7(t) < Cp[l + p] + o(l) 
as t — ► +oo. 

This theorem may be interpreted as a stability result for a perturbation of a homogeneous reaction- 
diffusion equation by the buoyancy coupling. The proof is based on the construction of super- and 
sub- solutions, and a bound on the decay of the solutions of advection-diffusion equations that is 
uniform in the advection flow. 

The third result of this paper deals with the Boussinesq system in a narrow domain. It has 
been shown in |14j that if a vertical strip is sufficiently narrow and gravity is sufficiently weak then 
solutions of the Cauchy data become planar as t — * +oo. The following theorem generalizes this 
result to inclined cylinders. 

Theorem 1.3 Let e = (ei, e<i) be the unit vector in the direction of gravity and let pj = pej, j = 1,2 
and let the initial data (To,uo) be as in Theorem \l. 6 A There exist two constants \q and po so that if 
the domain is sufficiently narrow: A < Ao and gravity is sufficiently small: p < po then the burning 
rate is bounded by 

V(t) < c + Cp 2 + o{\) as t +oo. (1.16) 
Moreover, the front is nearly planar in the sense that 

1 /"* 

N z (t) = - \\T z (s)\\ 2 2 ds < Cp\ + o(l) as t -» +oo. (1.17) 
t Jo 
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The main observation of this theorem is that only the gravity strength in the direction perpendicular 
to the strip enters in the upper bounds Q4.35|) and (|4.36l) . 

The paper is organized as follows: Theorem II .11 is proved in Sections [21 and 13 Theorems II .21 and 
11,31 are proved in Section 0J 
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2 The finite domain problem 

We consider in this section the approximating problem 

- cT x + u • VT = AT + f(T) (2.1) 
— cu x + u • Vu — erAu + Vp = pTe 

V • u = 0, 

in a finite domain D a = [—a,a] x x [0, X] y , a > 0, with the boundary conditions 

dT 

T(-a, z) = 1, T(a, z) = 0, — = at z = 0, A (2.2) 

and 

w = 0, u = at z = 0, A and v(±a, z) = u>(±a, z) = at x = dba. (2-3) 

One can show with the techniques of the present section that a solution T a , u a of (|2,lj) in D a with 
the boundary conditions (|2.2j) and (|2.3|) exists for all c£l. However, given an arbitrary c there is 
no way to control the limit of T a and u a as a — > oo. Hence, following the standard procedure, we 
impose an additional constraint (|l.llj) . This ensures that the non-trivial part of the solution does 
not escape to infinity when we pass to the limit a — > oo. 

Proposition 2.1 There exists a speed c a S R so that there exists a solution (T a ,u ) of \2.1]) in D a 
with the boundary conditions \2.°J\) and \2. S\) such that 

max T a (x,z) = e . (2.4) 

x>0,2G[0,A] 

We denote the corresponding solution as (c a ,T a ,u a ). Moreover, there exists ao > and a constant 
C > that is independent of a, so that we have for all a > oq 

\c a \ < C, (2.5) 

and 

/ \VT a \ 2 dxdz+ \Vu a \ 2 dxdz + Hujoo < C. (2.6) 

JD a JD a 

Moreover, the uniform Holder estimates hold: there exists ao > and a constant C > independent 
of a so that we have for all a > ao 

ll^alb^ODa) + II u a || CL" (A,) + H^oHci-"^) ^ C ( 2 - 7 ) 

provided that < a < 1. 
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Proof. The proof consists of two parts. First, we introduce a family of problems depending on a 
parameter r £ [0, 1] so that at r = we have a simple linear problem without advection or coupling 
and at r = 1 we have the full problem (|2.1j) with the correct boundary conditions. The normalization 
condition (|2,4j) is imposed for all r E [0,1]. We obtain the a priori bounds as in Q2.5J1 . Q2.6|) and 
(|2.7j) for such solutions that are uniform in r £ [0, 1]. In the second step we use the a priori bounds, 
the Leray-Schauder topological degree argument and the information on the linear problem at r = 
to show that solutions of the nonlinear coupled problem at r = 1 exist. We drop the subscript a 
throughout the proof to make the notation less cumbersome. 

Step 1. A priori bounds for solutions. Let us first define a one-parameter (homotopy) 
family of finite domain Boussinesq problems in the vorticity formulation 

-c T T T x + tu t ■ VT T = AT T + r/(T r ) (2.8) 
-c T u% + u r • Vw T - aAuj T = rpe • V X T := pr[e 2 T^ - e x T£\ 
u T = w r x ~ v T z , V • u r = 0. 

As mentioned above, r is the homotopy parameter: r £ [0, 1], with r = corresponding to the linear 
problem, and r = 1 to the full problem (|2.1[) - (|2.3|) . The problem (|2.8|) is posed in D a with the same 
boundary conditions 

= 0, w T = lo t = for z = 0, A (2.9) 

oz 

and 

T T (-a,z) = l, T T (a,z) = 0, v T (±a, z) = uj t = (±a, z) = for x = ±a, (2.10) 
as (|2.ip . We also require that 

maxT T (x,z) = 6 (2.11) 

x>0,z 

and obtain a priori bounds on c r , T T and uj t . We drop the superscript r below wherever it causes no 
confusion. The general plan is as follows. First, we bound the speed c above and below by a linear 
function of ||f ||oo i n Lemma 12.21 Next we bound HuH^ from above by a linear function of ||VT||2 
in Lemma 12.41 The other direction, a bound on ||VT||| in terms of a linear function of HuH^ is 
established in Lemmas 12.51 and 12 . 61 Since the latter bound is quadratic in ||VT||2, the last estimates 
allow to obtain a uniform bound on this quantity, from which all other a priori bounds follow in a 
fairly straightforward manner: see Corollary 12 . 71 and Lemma 12.81 

We begin with a lemma that bounds the speed c in terms of the horizontal flow velocity |H|z,°°(.D a ). 



Lemma 2.2 Let (c, T, u) satisfy \2.tfy - \2.1U\) with the normalization H2.ll)) and let u = (v, w). Then 
there exists ao > so that for all a > ao we have 

-1 - t|H|oo <c<1 + Mt + t|H|oo. (2.12) 

Proof. First, we observe that the function iPa(x) = Ae~ a<yX+a ^ is a super-solution for the reaction- 
diffusion-advection equation with the flow u fixed if A > 1 and 

Mr , 

c>a + +r \\v oo, (2.13) 

a 

that is, 

_A +TU .V^> AVu + t/(Va), (2-14) 
ox 
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provided that (|2.L-i|) holds with 

f(T) 

M= sup il_2. 

0<T<1 1 

Furthermore, we have 

T(-a,z) = l< A = ifj A (-a), T(o, z) = < if> A (a) (2.15) 

at the two ends of the domain D a . We now show that this together with (|2.14|) implies that 

T(x,z) <4> A (x) (2.16) 

for all (x, z) £ D a and A > 1. Indeed, consider the family of functions i/ja(x). Then all ipA are super- 
solutions in the sense that the inequality 1)2. 14|) holds. Moreover, as the maximum principle implies 
that < T < 1, for A > 5e 2aa sufficiently large we have tp A (x) > 5 > T(x, z) for all (x, z) £ D a . We 
define 

A = inf {A £ R : ip A (x) > T(x, z) for all (x, z) £ D a ] . 

The previous argument implies that Aq is finite, Aq < 5e 2aa and, moreover, clearly Aq > 0. Observe 
that since the domain D a is compact, we should have ipA {x) > T(x,z) - otherwise this inequality 
would be violated for A slightly larger than Aq at some point in D a . Moreover, the equation 
ipA Q (x) = T(x, z) should have a solution. We claim that Aq = 1. Indeed, otherwise the point (xo, ^o) 
that solves ipA (xo) = T(xq, Zq) cannot be at the boundary of D a because of the boundary conditions 
on the function T. Hence this point has to lie in the interior of D a . The continuity of i/ja(x) with 
respect to A implies that the graphs of tjjA (x) and T(x, z) are tangent at (xo, Zq). Then the strong 
maximum principle implies that i/ja ( x ) = T(x,z) which is a contradiction, as they differ on the 
boundary. Hence we conclude that Aq = 1 and thus Q2.16JI holds for all A > 1 and thus for A = 1, 
so that 

T{x,z) < e - a ( x+a \ (2.17) 

However, the existence of such a super-solution contradicts the normalization condition 1)2.11(1 
if a > ln(0Q 1 )/a because (|2.11JI implies that there exists zq so that T(0,Zq) = 9q. Therefore, the 
existence of a solution T that satisfies (|2.1ip implies 

/ Mt\ 

c< inf aH + t\\v\\oo < 1 + Mr + t|H|oo (2.18) 

a>ln(0o x )/a V a J 

provided that a > ln(l/^o)- This proves the upper bound in (|2.12|) . In order to prove the lower 
bound we observe that the function <$> = 1 — ^ al \ x ~ a ) is a sub-solution for T with the flow u fixed if 

c < -a - tWvWoo. (2.19) 

That is, if (12T91) holds, then T(x,z) > 1 - e a ( x - a \ This is shown in a way similar to the proof of 
(|2.17|) under the assumption (|2.L-{|) above. However, </>(0) = 1 — e~ aa > 9$ for 

ln((l - ^o)" 1 ) 

a > — ^ J — L . (2.20) 

a 

This implies that max x >o T(x, z) > 0(0) > 8q provided that both (|2.19|) and (|2.20|) hold. Hence in 
order for (|2.11|) to be possible we need 

c> sup [— a — rH^Hoo] > — 1 — t||u[|oo (2-21) 
a> M(irMzl) 
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provided that a > ln((l — ^o)" 1 )- This is the lower bound in (j2.12|) and the proof of Lemma 12.21 is 
complete. □ 

Next we establish a bound on ||u||2/x>(£) o ) and H^H^oo^) in terms of || S/T\\ L 2r Da y These bounds 
are all obtained from the following type of estimates. 

Lemma 2.3 Let S a = [—a, a] x xfl y be a finite cylinder with a smooth bounded cross-section Q G M. d , 
d = 1,2. Let <j> be a function that satisfies either of the following three conditions: (i) 4>(x,y) = 
on the whole boundary dS a , (ii) <j)(x,y) = for y € dQ and 9 ^Q x ' y ^ = for x = —a, a, or (Hi) 

gin = f or V £ <9$7, and (j)(x, y) = for x = —a, a. Then there exists a constant C that depends 
only on the domain £1, but not on the cylinder length a, so that we have 



(Sa) <C[||A^|| L2(5a) + |H| i2(Sa) ]. (2.22) 



Proof. Let Q be any cylinder of the form [xq, xq + 1]x£1 C S a with —a < xq < a— 1. The standard 
interior elliptic estimates up to the boundary [20] can be applied to Q in all the three cases (i)-(iii). 
The corners at x = ±a are not an obstacle. Indeed, both in the case of the Dirichlet and Neumann 
boundary conditions prescribed on the lines x = ±a, one can extend the solution to a larger cylinder 
[—a — 1, a + 1] X O by reflecting the solution across the line x = ±a, either in the even or odd 
way, respectively. Hence the usual elliptic estimates up to the boundary can be applied to all such 
cylinders Q to obtain 

U\\ miQ) < C [IIA^I^j + U\\ L 2 {Sa) ] (2.23) 
in all three cases (i)-(iii). Then the Sobolev embedding theorem in dimensions d = 2, 3 implies that 

U\\l°°(q) < CU\\ H 2 {Q) < C [\\A<t>\\ L2{Q) + U\\ L 2 (Q) ] < C [||A^|| l2(5q) + U\\ L 2 (Sa) ] 

with the constant C that depends only on the domain Q. □ 

This lemma can be easily extended to higher dimensions using the appropriate Sobolev embed- 
dings. It implies immediately the following bounds on ||u||oo and ||a>||oo in terms of || VT\\ L 2^ Da y 

Lemma 2.4 Let (c, T, u) satisfy \2. 8\) - i2.1(J\) with the normalization \2.11\) . There exists oq > 
and a constant C > so that, for all a > oq, it holds that 

[|u|U~ (Da) < C\\VT\\ L 2 (Da) . (2.24) 

and 

IMU» (Da) < C||VT|| L2(Da) [1 + ||Vr|| i2(Ai) ] . (2.25) 
Moreover, Vu satisfies the same bound: 

||Vu|| ioo(Da) < C\\VT\\ L 2 {Da) [1 + ||Vr|| i2(Da) ] . (2.26) 

Proof. We use the vorticity equation 

-cuj x + u • Vw - aAu = pr(e • V^T), uj = on dD a . (2.27) 

Case (i) of Lemma 12.31 implies that 

IMU-tA,) < C [\\VT\\ L 2 {Da) + (|c| + ||u||oo)||Va;|| jL 2 (jDa) + |M| L2{Da) ] . (2.28) 
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Here the constant C depends only on p and A. Note that multiplying the vorticity equation by uj 
and integrating by parts, using the boundary conditions we obtain 

\X7uj\ 2 dxdz = rp J • V^T^ udxdz < Tp|| VT||2||w||2- 

The Dirichlet boundary conditions for uj imply that the Poincare inequality applies to uj so that 
IMlL 2 (D a ) < (^/ 7r )ll^ w llL 2 (D a )- Hence we obtain 

||Vw|| 2 < -rp||Vr|| 2 (2.29) 

and thus 

Ml^Do) < C||VT|| i2(Da) , (2.30) 

with the constant C independent of the cylinder length a. This, together with (|2.28|) and the bound 
(|2. 12|) on the speed c implies (|2.25|) . provided that we show ()2.24j) . 

We now prove (|2.24|) . The horizontal flow component v satisfies the Poisson equation 

dv 

Av = -uj z , v(±a,z) = 0, ^ = 0, atz = 0,A. (2.31) 

The boundary conditions at z = 0, A are obtained from v z = w x — uj = as follows from (|2.9|h The 
third case (hi) of Lemma 12 . 31 implies that 

||«IU<»(D a ) < C[\\SJu\\ L 2 {Da) + \\v\\ L 2 {Da) }. (2.32) 

The first term in the right side is bounded by (|2.29() . In order to bound the second we multiply 
(|2.31j) by v and integrate to obtain, using the boundary conditions and (|2.29|) 

Vv\ 2 dxdz = I w z {x,z)v{x,z)dxdz <\\uo z \\2\\v\\2 <C\\\/T\\2\\v\\2- (2.33) 

D a J D a 

Now, observe that (|2.31|) . the Neumann boundary conditions for v and the Dirichlet boundary 
condition for to at z = 0, A imply that 



d 2 f 

/ v(x, z)dz = 0. 
dx z J 



It follows then from the Dirichlet boundary conditions for v at x = ±a that 



v(x,z)dz = (2.34) 



o 



for all x. One may alternatively deduce Q2.34|) from incompressibility of the flow u and the boundary 
conditions. Therefore, it follows from the Poincare inequality that |M|L 2 (D a ) < (A/27r)||Vt>|| i 2(£ )a ). 
Thus, ()2.33j) implies that both ||Vf Hz^^) < C\\ VT\\ L 2( Da j and ||«||£2m ) < C\\ VT\\ L 2^ Da ^ with a 
constant independent of a. Hence, ()2.32(l implies (j2.24|) for the horizontal flow component. 
The vertical flow component satisfies 

dw 

Aw = lu x , w(x, 0) = w(x, A) = 0, (±a,z) = 0. (2.35) 

The Neumann boundary condition at x = ±o is deduced from the relation w x = u + v z and the 
Dirichlet boundary conditions for v and ui at x = ±a. The case (ii) in Lemma 12 . 31 implies that 



\w 



\L^{Da) < C[\\Vu\\ L 2(D a ) + \H\LHDa)}- ( 2 - 36 ) 
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As before, we use (|2.29j) to bound the first term in the right side. In order to bound the second we 
multiply (|2.35f) by w and integrate, using the boundary conditions and (|2.29j) again, to obtain that 

/ | Vw\ 2 dxdz = — u> x (x, z)w(x, z)dxdz < \\uj x \\2\\w\\2 < C\\ VT||2||w||2- (2.37) 

J D a J D a 

The Dirichlet boundary conditions for w at z = 0, A imply that ||^||i, 2 (D a ) — / V 7r ll^ 7u, llL 2 (D £l )- Thus, 
(|2.37f) implies that 

\\Vw\\ L 2 {Da) < C\\VT\\v (Da) , (2.38) 

and hence ||iw||£2m) < C\\ VT\\ L 2^ Da ^ with a constant independent of a. Therefore, now (|2.36f) 
implies (|2.24j) for the vertical flow component. Thus, the proof of ()2.24j) is complete. We recall that 
then (|2.25|) follows as well, as explained in the paragraph below (|2.3()j) . 

In order to complete the proof of Lemma 12.41 it remains to bound the derivatives of u. First, we 
observe that the function ip = v z satisfies the boundary value problem 

-Alp = uj zz , V = on dD a . (2.39) 

Hence, case (i) of Lemma 12.31 applies to the function ip. Moreover, the elliptic estimates for u, as 
in (|2.23j) imply that ||^|| L 2( Da) < || Au; || L 2( D a) < C||Vr|| L 2 (Da )(l + ||Vr|| L 2( Da) ). Hence, the same 
proof as in the derivation of the bound (|2.25j) applies to ip and we obtain that 

K|Uoo ( A0 < C||VT|| i2(Da) (l + ||VT|| L 2 (Dq) ). 

This, together with ()2.25j) implies that 

IK|| L oc (Da) < c||vr|| L 2 (Da) (i + ||vr|| i2(D[!) ). 

The other pair of derivatives, v x and w z , do not satisfy a homogeneous boundary condition on the 
lines x = ±a. Therefore, one cannot apply the standard elliptic estimates up to the boundary to the 
function ry = w z = —v x (the second equality follows from the incompressibility of the flow). In order 
to circumvent this difficulty, we extend the function w to a larger cylinder D a+ \ = [— a— 1, a+1] x [0, A] 
by setting w(—a — x,z) = w(—a + x, z) and w(a + x, z) = w(a — x, z) for < x < 1. The resulting 
function is of a class C 2 (D a+ i) since w(x, z) satisfies the Neumann boundary condition at x = ±a. 
This also extends the function r] = w z to the larger cylinder. Moreover, r/ satisfies the Neumann 
boundary condition along the horizontal lines z = 0, A: 

Vz = Wzz = ~v zx = on z = 0, A, 

and 

At? = u xz , 

with the function ui extended to the larger cylinder by the same reflection. Hence, the interior elliptic 
estimates up to the boundary for solutions of the Neumann problem imply that 

IMIiP ( Q) < \\Ar]\\ L 2 {Da) + \\r]\\ L 2 {Da) 

for any rectangle Q = [xq, xq + 1] x [0, A] that is strictly contained inside the larger cylinder D a+ \. 
Therefore, the Sobolev embedding theorem together with the above estimates imply that 

h\\ L ~>(D a ) < C[\\Ari\\ L 2 (Da) + \\r]\\ L 2 (Da) ] = C[\\uj xz \\ L 2 {Da) + |M| L 2 (Da) ]. (2.40) 
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However, as the function uj satisfies the Dirichlet boundary conditions in D a , we can apply the 
estimate ()2.23j) to the function u up to boundary, to obtain 



UJ 



\jp(D a ) < C [\\ AuJ \\mD a ) + IMU 2 (D a )] • 



We now use the vorticity equation (|2,27|) to bound ||Aw||2^(d ) an d the estimate (|2.30j) to estimate 
ll w lll, 2 (D a )> and conclude that 

\\u xz \\ L 2 {Da ) < C||VT|| i2(Da) (l + ||VT|| L2(Da) ). (2.41) 

Furthermore, the estimate (|2.38|) for ||Vio||x,2/£) ) implies that 

h\\L*(D a ) = \\w z \\mD a ) < l|v^|| L 2 (Da) < c||vr|| L2(£ , a) . (2.42) 



We infer from the bounds (j2.40[) . (|2.41jl and (|2.42j) that 

NIU»(£»„) < C||Vr|| i2(Da) (l + \\VT\\ L2{Da) ). 



This proves the uniform bound on w z and hence the proof of Lemma 12.41 is complete. □ 

Let us now proceed to estimate ||VT[|^2m o ) in terms of \\v \\L°°(D a )> a bound in the direction 
opposite to that in Lemma 12.41 Most importantly, we will bound the square ||Vr||| in terms of a 
linear function of II 

t^||oo* As we are unable to obtain such bound by the standard ellitpic estimates, 
we have to proceed with an explicit calculation. As a preliminary step we show the following. 

Lemma 2.5 Let (c, T, u) satisfy \2.£\ )- \2.1(J\) with the normalization \2.11\) . Then there exists a 
constant C > and a constant oq > so that we have for all a > oq and < r < 1 

I \VT\ 2 dxdz+ [ T x (a, z)dz < C [1 + \\vWoo] . (2.43) 

JD a JO 

Proof. Recall that the function T satisfies 

—cT x + ru • VT = AT + r/(T) (2.44) 

with the boundary conditions 

dT 

T(-a, z) = 1, T(a, z) = 0, — = at z = 0, A. (2.45) 

We multiply (|2.44() by (1 — T) and use the boundary conditions and incompressibility of the flow u 
to obtain 

cA 



-= I T x 



f T x (a,z)dz+ [ \VT\ 2 dxdz + r f (I - T)f{T)dxdz. (2.46) 
Jo Jo a J 

Hence Lemma 12.21 and the fact that (1 — T)f(T) > imply that 

/ \VT\ 2 dxdz+ [ T x (a, z)dz < ^ < C [1 + |M|oc] (2.47) 

JD a JO 2 

and Lemma 12.51 is proved. □ 

In order to close the bounds ()2.12|) . ()2.24|) and (|2.43l) we need to bound the integral of T x in 
(|2.43[) . This is done in the next Lemma. 
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Lemma 2.6 Let (c, T, u) satisfy \2. with the normalization h2.11\) . There exists a constant 
C > and a constant oq so that we have for all a > ao and < r < 1 

0<- / T x (a, z)dz < C [1 + \\VT\\ 2 ] . (2.48) 
Jo 

f X 

Proof. In order to find a bound for / T x (x = ±a, z)dz we introduce 

J o 

1 f x 

I( x ) = t / T(x,z)dz 
* Jo 

and integrate equation 1|2.8[) for T m z. Using the boundary conditions we obtain 

-L xx = G(x), J(-o) = 1, /(a) = 0, G(x) = jj f(T(x, z))dz - J (tu ■ VT - cT,)y . (2.49) 
This equation can be solved explicitly: 

/X 
{x - s)G{s)ds + Ax + B 
-a 

with constants 

A = ~Ya + YaF {a ~ s ^ ds > B = \ + \f i ' a - s ^ ds 

•J CL o — Q, 

that are determined from the boundary conditions. Thus, we have 

I x (-a) = A, I x (a) =A-f G(s)ds. 

J —a 

Using the expression for the function G{x) in (|2.49[) . we now infer that 

< -4(a) = Y a + Ya f a {a + s)G{s)dS = h + hf a / A( ° + X)/(T(X ' Z)) ^T 
t f a [ x , , „ m/ ,dzdx c f a f x , ^ dzdx 

Integrating by parts, using the boundary conditions and incompressibility of u we obtain 
r f a f x , ,dzdx c f a f x dzdx 



+ ^Lh v(x '* )T(x ' z) —-^Lh T a ■ 

However, the normalization condition (|2. implies that f(T(x,z)) = for x > since there is no 
reaction to the right of x = 0. Therefore, we can drop the third term above. This is one of the 
crucial points in the proof of the current lemma. Hence, we conclude that 



v(x, z)T(x, z)- 



c f a f x rT1 dzdx 1 t f a f x rtriU ..dxdx , 

WJ. T —^ + 2hL /m.,.))— + r|MU + W . (2.50) 
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We used the fact that < T < 1 to bound the last term above. Next, we look at I x {—a): 



< -I x (-a) 



1 1 

2a ~ 2a 



(a — s)G(s)ds 



1 r 
2a ~ 2a 



-a JO 



x dzdx 
(a - x)f(T(x,z))—^- 



— / / (a — x)u ■ VT(x, z)dzdx / / 

2a J - a Jo 2a J_ a J 



a />A 



(a - x)^ 



dzdx 



We can drop the second term above, as (a — x)f(T) > 0, so that, after integration by parts, we get 



< -4(-a) < — + — / / t;x,z)Tx,z) — 

1 

< 7T + r M U + C • 
2a 

Let us now put together (|23U)) . (|23T|) and ((2~1H)) We observe that, with F = r f /(T) 
have the following three inequalities: 

c = I x {a)-I x {-a)+F 
1 F 

< -I x {a) < — + — + t||v||oo + |c| 
2a 2 

< -Ix(-a) < — + t||u||oo + \c\- 



a f x rT1 dzdx , 
, / T— r — -t- | c| 
2a ./-a Jo A 



(2.51) 

dxdz 



A ' 



wc 



This implies that 



and thus 



Lemma 12 . 21 implies then that 



F < - + 2t[|'u[| oq + 4|c| 
a 



-4(a) < 3|c| H h 2r||f Hoc for a > ao 



(2.52) 



-4(a) < C 

Finally, Lemma 12.41 implies that 

-4(a) < C 



i + tIMIoo + - 



i+ vr 2 + - 



for a > ao- 



Thus, Lemma 12.61 is proved. □ 

The previous lemmas imply uniform bounds that we summarize as follows. 

Corollary 2.7 Let (c, T, u) satisfy \2.tfy - \2.10}) with the normalization \2.11\) . There exists a con- 
stant C > and ao > so that we have for all a > ao and < t < 1 

ll u llL°°(D a ) + l|Vu|| L oo (Da) + \\uj\\ L 2 {Da) + |c| + || Vr|| L 2 (Da) +T f(T)dxdz < C. (2.53) 

JD a 

In particular, as a consequence we also have 



I Vu IIl2(d ) + IMIl^a,) < c. 



(2.54) 
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Proof. Lemmas 12.61 and 12.51 imply that 

J \VT\ 2 dxdz < C [1 + rlMloc + ||VT|| L 2 (Da) ] 



Then Lemma 12 . 41 implies that 



||vr||i 2(I3o) <c(i + ||vr|| L2(Da) ) 



and thus the estimate on || VT|| i 2(£> a ) in (|2,53|) holds. Then Lemma 12.41 implies the bounds on 
ll u llz,°°(D a )) ll^ u llL°°(D a ) an d ||^||l2(d )- The bound on |c| in (|2.53j) now follows from Lemma 12.21 
Finally, the estimate on the total reaction rate follows from the above bounds and (|2.52|) . One can 
elliminate the factor r in front of the total reaction rate in (|2.53[) : actually, one can show that it 
remains bounded as r — > 0. However, unlike the other estimates in (|2.53j) . we will use the bound on 
the total reaction rate only at r = 1. □ 

It remains to prove the uniform Holder (^'"-estimates for T(x, z), to and u in order to finish the 
proof of Proposition 12.11 

Lemma 2.8 There exist two constants C > and ao > so that the following bound holds for all 
a > oq: 

IMIcv^A.) + llullc^CDa) + ll r llc ,1 ' 0! (£»a) - C (2.55) 

provided that < a < 1. 

Proof. The bound for T follows from the standard elliptic local regularity estimates up to the 
boundary [20], the C 1 -bound on the flow u and the uniform bound on the speed c in Corollary 12.71 
The Holder estimate for u follows then from the vorticity equation (|2.27|) with the Dirichlet boundary 
conditions, the above mentioned C 1,Q -bound on T, the same uniform estimates in Corollary 12 . 71 and 
the same results of |2Uj . Finally, the Holder bounds on u follow from the Poisson equations (|2.31f) 
and (|2,35|) on the horizontal and vertical flow components, respectively, and the Holder estimates 
for uj. □ 

This completes the proof of the a priori bounds in Proposition 12. II We now turn to the proof of 
the existence part of this proposition. 

Step 2. The degree argument. The a priori bounds proved in the first step of the proof 
allow us to use the Leray-Schauder topological degree argument to establish existence of solutions 
to the problem H2.8|) - ()2.10|) with the normalization (|2.11|) in the bounded domain D a . This method 
of construction of traveling wave solutions goes back to We introduce a map 

K T : (c,oj,T) -► (8 T ,n T ,Z T ) 

as the solution operator of the linear system 

- cZl + ru • VZ T = AZ T + r/(T) (2.56) 
-cn T x + u • VST - aAQ T = T P [e 2 T x - e x T z ] 



in D a with the no stress boundary conditions 

dZ 1 



dz 

and 



0, w T = n T = at z = 0, A (2.57) 



Z T (-a,y) = l, Z T (a,y) = 0, u T = n T = at x = ±a. (2.58) 
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Here the unknown flow u r = (u T , w r ) and the given flow u are the incompressible flows corresponding 
to the vorticities ft T and uj, respectively, and satisfying the no-stress boundary conditions. The 
number 9 T is defined by 

8 T = 6*n — maxTfx, z) + c. 

The operator K, T is a mapping of the Banach space I = Kx C 1,a {D a ) x C l,a {D a ), equipped with 
the norm ||(c, uj,T)\\x = max(|c|, H^H^am^, ||r||cn,am o ^), onto itself. A solution q T = (c T , u> T , T T ) 
of H2.8|) - (|2.1U|) is a fixed point of IC T and satisfies /C T q r = q r , and vice versa: a fixed point of K, T 
provides a solution to (|2.8j) - ([2.10|) . Hence, in order to show that (j2.8j) - (|2.10|) has a traveling front 
solution it suffices to show that the kernel of the operator T T = Id — K, T is not trivial. The standard 
elliptic regularity results in [2U| imply that the operator JC T is compact and depends continuously 
on the parameter r G [0, 1]. Thus the Leray-Schauder topological degree theory can be applied. Let 
us introduce a ball Bm = {||(c, u>, T)\\x < M}. Then Lemma 12.81 and Lemma 12.21 show that the 
operator T T does not vanish on the boundary 8Bm with M sufficiently large for any r £ [0, 1]. It 
remains only to show that the degree deg(^i, Bm, 0) in Bm is not zero. However, the homotopy 
invariance property of the degree implies that deg(^v, Bm,0) = deg(.Fo, %, 0) for all r £ [0,1]. 
Moreover, the degree at r = can be computed explicitly as the operator Tq is given by 



Jb(c,w,T) = (maxT(x,y) - 9 ,uj,T - T c ; 

x>0 



Here the function T${x) solves 



dx 2 dx 

and is given by 



°+c^ = 0, T c (-a) = l, T » = 



ca p—ca 



The mapping T§ is homotopic to 



$(c,w,T) = (maxT o c (x,y)-0 o ,w,T-T o c ) 

that in turn is homotopic to 

§(c,w,T) = (T o c (O)-0 o ,^T-T o cS ), 

where is the unique number so that Tq*(0) = #o- The degree of the mapping 6 is the product of 
the degrees of each component. The last two have degree equal to one, and the first to —1, as the 
function Tq(0) is decreasing in c. Thus deg^o = — 1 and hence deg.Fi = — 1 so that the kernel of 
Id — /Ci is not empty. This finishes the proof of Proposition 12.11 □ 

Remark 2.9 Observe that the C 1 ' -regularity ofT, u and to can be bootstrapped to C 2 ' a -regularity: 
we have 

IMIc 2 >«(D a ) + ll u llc 2 . Q (D a ) + ll^llc^A,) < C ( 2 - 59 ) 

provided that < a < 1 . 
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3 Identification of the limit 



In order to finish the proof of Theorem 11.11 we consider the solutions (c a ,T a ,u a ) constructed in 
Pr op osition 12 . 1 1 and pass to the limit a — ► +oo. The a priori estimates in the same proposition imply 
that we can choose a subsequence a n — > oo so that T n (x,z) = T an (x,z) converges uniformly on 
compact sets to a function T(x, z), while the flow u n (x, z) = u an (x, z) converges to a flow u(x, z) = 
(v, w) and the front speeds also converge: c n = c an — ► c. The vorticity functions co n (x, z) = uj an (x, z) 
converge to the limit uj = w x — v z . The limits satisfy the uniform bounds 



c + u 



+ IHU + J \VT\ 2 dxdz + j f(T)dxdz + J \Vu\ 2 dxdz < C (3.1) 



that follow from Corollary 12.71 and the Holder estimates Q2.7|) and Q2.59JI . The regularity estimates 
on (T a ,u a ) imply that the limit functions T and u satisfy the Boussinesq system 

-cT x + u • VT = AT + f(T) 

-coj x + u-Vuj- aAuj = p(e ■ V^T) 

oj = w x -v z . (3.2) 

Moreover, the boundary conditions on the lateral boundaries hold for T and u: 

dT 

-7— = 0, w = uj = on z = 0, A. (3.3) 

oz 



The normalization condition 



maxT(x,z) = 6 (3.4) 

x>0 



is also satisfied. 

Therefore, to finish the proof of Theorem II. H it remains only to show that (i) T converges to 
a constant 9- as x — * —00 and T — > as x — > +00, (ii) u — * as x — > +00, and (iii) 9— = 1 
if the reaction rate satisfies f(T) < (T — ^o)+/A 2 . First, we note that the uniform L 2 -bound on 
VT in ()3.1|) implies that T converges to two constants 9- and + as x — > ±00, possibly passing to 
a subsequence x n — > ±00. The elliptic regularity results imply that actually T converges to these 

— - ■ ^ - — - - - - - - — - / ™ - » — 

that f(9-) = f(9+) = 0. Furthermore, integrating (|3.2j) we obtain 

c(0_-0 + ) = Jf(T)^. (3.5) 

In order to identify the limits 9± we will make use of the following lemmas that provide some 
additional information on solutions on a finite domain before the passage to the limit. The first 
result describes the behavior near the right end x = a n . 



Lemma 3.1 There exists a sequence a n — > 00 so that 

dT n (a n ,z) 



Ox 

as n — > 00, uniformly in z. Moreover, we have lim n ^ 00 T n (a n — xq, z) =0 for all xq 6 



(3.6) 
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Proof. We introduce a shifted solution <£ n (x, z) = T n (x + a n , z) , v n = u n (x + a n , z) denned in the 
domain — 2a n < x < 0. The functions $ n and v n satisfy the same a priori bounds Q3.1JI as T n and 
u n and hence they converge as n — * oo to some limits $ and v that satisfy 

-c& x + v • V$ = A$, $(0,z) = 0, x<0 (3.7) 

as f(& n ) = for x > — a n and thus in the limit /(3>) = 0. The function <I> satisfies the Neumann 
boundary conditions at z = 0, A. The uniform upper bound on ||V$ n ||2 together with the elliptic 
regularity results imply that $ has to converge to a constant 3>_ as x — ► — oo along a subsequence. 
We note that, as < <J>(x, z) < 6>o, the constant 3>_ satisfies the same bounds: 

< < 6» . 

Integrating 1)3. 7 Jl we obtain 

cA$_ = J<f> x (0,z)dz < 0. (3.8) 

Hence, either 3>_ = or c < 0. In the former case <I> = and hence $ x (0, z) = for all z. That 
implies that both T™(a n ,z) — > as n — ► oo and T™(a n — £o,z) = c ^n(— ^O;- 2 ) -> as n -> oo, as 
claimed in Lemma 13. II It remains to rule out the second case, c < 0. This is done in the next lemma 
that provides a crucial lower bound on the speed c n . In particular it shows that c > - this will 
conclude the proof of Lemma 13.11 

Lemma 3.2 The front speed is positive, c > 0. 

Proof. Integrating the temperature equation in (|3.2|) for T n , we obtain 

°n^= [ ^ L {a n ,z)dz - [ ^ !1 (-a n ,z)dz+ [ f(T n )dxdz> [ ^ L (a n ,z)dz + f f(T n )dxdz. 



ax Jo ^ x J Jo 

(3.9) 

Observe also that we have a uniform bound 

J f(T n )dxdz\ (^J \VT n \ 2 dxdz^j > C (3.10) 

that follows from the fact that T n (0,z) < 9q an d T n (—a n ,z) = 1. The proof of 1)3.10)) is as in |13j : 
there exists zq G (0, A) such that both 

o 

\VT n (x,z )\ 2 dx < 3 / |VT n | 2 ^— , D n = [-a n ,a n ] x x [0, A] z , 

and 

o 

f(T n (x,zo))dx < 3 / /(T n (x,z)) — — . 



— an -Dn 

Let xi be the left-most point so that T n (x\, Zq) = 1 — : 

f 1 — 6 

Xi = inf (if (-On, 0) : T n (x, z ) = 1 
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^ Q 

and x 2 > x\ be the left-most point so that T n (x2, Zq) = 9 H — : 

1 



X2 = inf < x G (-a„, 0) : T n (x, z Q ) = O + 



Existence of xi and x 2 is guaranteed by the fact that T n (—a n ,z) = 1 and T n (0,z) < 0o f° r a U 
z G [0, A]. Then the reaction rate f(T n (x, zq)) > C for x\ < x < x 2 so that 



/•X2 r 

C\xi-x 2 \< / /(T n (a;,zo))^ < 3 / f(T n (x,z)) x 

and 

(l-0o) 2 f°*,„ m , m2 - f ^ adxdz 



4|xi — x 2 | 



< / 2 |V^n(x,z )| 2 dx<3 / |VT n | 



Multiplying these two inequalities, we arrive at (j3,l()|) . 

The estimate (|3,10j) and the uniform upper bound on ||VT n || 2 in Corollary 12 , 71 imply that 

/ f(T n )dxdz > C. (3.11) 

JDn 

Then, passing to the limit in (|3.9|) . and using (|3.8|) we obtain 

cA(l - > C > 0, 



as 



(a n ,z)dz -> / ^(O^ete, 



9x 

with the function $ as in the proof of Lemma 13.11 Now, we recall that <£_ < 8q < 1 and thus the 
front speed c > 0. This finishes the proof of Lemma 13.21 and hence also that of Lemma 13. II □ 

Lemma 13.21 and ()3.5|) imply that 0_ > 9 + . However, if 0_ = + we have f(T) = everywhere 
and hence (|3,2j) is a linear equation. The maximum principle implies that T = const in this case. 
The last condition in (|3.1j) implies that this constant has to be equal to 0o- Hence, either 9- > 9 + 
or T = O . 

Let us now rule out the special case that 0_ = 9 + = 9q. 
Lemma 3.3 The left and right limits 9- and 9+ satisfy 9- > + . 
Proof. We have already shown that 9- > 0+ and, moreover, if 0_ = + then 

9_ = 9 + = O . (3.12) 

Hence, it suffices to show that the latter is impossible. Let us assume that ()3.12|) holds. As we have 
explained above, then 

T n — ► 0o, and dT n /dx — > uniformly on compact sets. (3.13) 
Then, integrating the equation 

dT 

-Cn-g^- + u n ■ VT n = AT n + f(T n ) 
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between x = and x = a n we obtain, as f(T n ) = in this region, 

f x f x f x dT f x dT 

c n T n (0,z)dz- v n (0,z)T n (0,z)dz= —^{a n ,z)dz- —^(0,z)dz. (3.14) 



JO 



o o 1 ^ Jo 



We now pass to the limit n — > oo in (|3.14|) . The first term on the left converges to c9q\, as we 
have assumed that T converges uniformly to 9q on compact intervals. The second term on the left 
converges to 

A 

v(0,z)e dz = 0, 



o 



as incompressibility of the flow u n and the boundary conditions at x = ±a imply that 

v n (0, z)dz = 0. 



The limit ()3.6j) in Lemma 13 1 1 1 implies that the first term on the right side of (|3.14|) converges to zero. 
Finally, the last term on the right side of (|3.14|) converges to zero because of (|3.13j) . Therefore, we 
obtain 

c\9 = 0. 

However, this implies that c = which contradicts Lemma 13.21 Hence, the case 9- = Q + = 9q is 
ruled out and thus 9- > 9 + . □ 

We continue the analysis of the behavior of the solution at the right end of the domain. 

Lemma 3.4 The gradient VT n converges to zero "as x — * +oo " uniformly in n, that is, for every 
e > there exists N £ N and R so that |VT ra (x, z)\ < e for all n > N and all R < x < a n . 

Proof. Let us assume that this is not the case. Then there exists £o > and a sequence b n — > +oo 
so that |VT n (6 n , z n )\ > eo for some z n G [0, A]. Note that Lemma 13. II implies that 

\a n ~ K\ oo- (3.15) 

Let us define the shifted solution *& n (x, z) = T n (x — b n , z), v n (x, z) = u n (x — b n , z) on the domain 
x € [—a n — b n , a n — b n \. Then ^ n and v n satisfy the same uniform bounds as T n and u n and thus they 
converge to a pair of functions ^f, v uniformly on compact intervals, together with their derivatives. 
The functions ^ and v are defined on the whole real line because of (|3.15|) . Moreover, the function 
^ has left and right limits *$>± as x — > ±oo. Hence, the same argument as in the proof of Lemma 
13. II shows that \I / must be equal a constant, as it has left and right limits and satisfies 

-c^ x + v • = A^. 

However, this contradicts the fact that max z |VW(0, 2r)| > eo- d 

The decay of the gradient of T n implies that the flow ahead of the front goes to zero for large x, 
uniformly in n. 

Lemma 3.5 The flow u(x, z) converges to zero on the right uniformly in n, that is, for any e > 
there exists R > and N £ N so that \u(x, z)\ < e for all R < x < a n . 

Proof. We choose N and R so that |VT n | < e for all n > N and x £ [R,a n ]. Then, we decompose 
T n = T l n n + T° ut with supp T l n n C {x < R + 1} and supp T° ut c{x>R}. We also require that both 
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T^ n and T° ut satisfy the same uniform gradient bounds as T n . Moreover, we have |VT° U *| < e. We 
also split uj n = lo™ + w™' and u n = u^ 1 + u° u * accordingly: 

/ dT in 1)1'"' \ 
-c n uC + u n • V<* - aAuC = P i e 2 ~ e i J > <C = an d£> an , 

and similarly for w™ 4 . 

We now bound |u^ n | and |u° u *| separately for sufficiently large x. First, we look at u^ n = 
(v™,w™). The function u;^ n satisfies a homogeneous equation 

~ Cn ~dx L ~ + Un ' ^ " = ° (3 - 16) 

in the rectangle Dn + 2 An = {R + 2 < x < a n , < z < A}, as T^ n vanishes in Dn + 2 t a„- The function 
satisfies a uniform C 2,a -bound - this is shown in the same way as the C 2,Q -bound for the full 
vorticity function oj in (|2.59|) . This in turn implies that the function ip(z) = uj™'(R+2, z) is uniformly 
bounded in C 2 [0, A]. Let g(x) be a smooth monotonic and positive cut-off function so that 

g(x) = 1 for R + 2 < x < R + 3 and g(x) = for x > R + 4. (3.17) 

Then the function can be decomposed as 

u™(x,z) = ^(z)g(x)+( n . 

The function £ n satisfies 

~° n ~dx + Un ' ~~ = ^" in DR+2 ' an > = on dD R+2,°- n - (3.18) 

The right side / n is given by 

f n := aip"(z)g(x) + aip(z)g"(x) - c n ?p(z)g'(x) - v n ip(z)g'(x) - w n ip' (z)g(x). 

It is supported ini? + 2<x<i? + 4 and is uniformly bounded since [|w^ l [|c 2 >«(D ) — C. Let us 
choose a > sufficiently small, then the function £ n (x, z) = ( n (x,z)e ax satisfies 

~ Cn ~dx + a °^ n + Un ' V ^ n ~~ aVn ^ n ~ aA ^ n + 2aa ~dx ~ aa2 ^ n = 9n in Dr + 2 <°™ > £n = on dD R+2 ,a n 

(3.19) 

with g n = f n (x)e ax . Multiplying (|3.19[) by £ n and integrating by parts, using the boundary condi- 
tions, we obtain 

a \^£,n\ 2 dxdz + (ca - a\\v\\oo - aa 2 ) / \£ n \ 2 dxdz < \\g n hUnh- (3.20) 

^ D R+2,a n JD R+2t a n 

However, as the function £ n vanishes at z = 0, A, the Poincare inequality implies that 

■2 



/ \V£ n \ 2 dxdz > ^ / \£ n \ 2 dxdz. 



J R+2,a n JJJ R+2,an 

Hence, the following upper bound holds 



S D 



\^ n \ 2 dxdz < \\g n \\l < C, 



R+2,a„ 
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provided that a is sufficiently small, since H^Hoo — C . Using (|3.2()|) onc6 again we conclude that 

\V£ n \ 2 dxdz < C. 



L 



DR+2,a r , 

Therefore, the function ( n satisfies 



/ [|VCn| 2 + |Cn| 2 ] e 2aX dxdz < C. 



This, in turn implies the same bound for the function 

ra n /"A i-a n r\ 

/ / \u™\ 2 e 2ax dxdz + / / \Vu™\ 2 e 2ax dxdz < C. (3.21) 
Jr+2 Jo Jr+2 Jo 

It follows that the L 2 -norm of decays uniformly in n: 

/ / \uj™\ 2 dxdz < e~ ar ° / / \uj™\ 2 e 2ax dxdz <Ce- ar °. (3.22) 

Jr Jo Jr Jo 

for vq > R + 5, and the same bound holds for Va->™: 



ra. n rX ra„ r-\ 

/ / \Vu™\ 2 dxdz < e~ ar ° / / \V uj™\ 2 e 2ax dxdz < Ce~ aT °. (3.23) 

Jr Jo Jr Jo 



As the function of™ satisfies the homogeneous equation (|3.16|) for x > R + 1 with a bounded flow 
u, the standard local elliptic estimates now imply that 

\ui™(x, z)\ < Ce~ ax for x > R + 5. (3.24) 

The W 2,p elliptic estimates imply then the uniform decay of the gradient of to™: 

\Vu™{x, z)\ < Ce~ ax for x > R + 5. (3.25) 

Now we can bound the flow = (v™-,w™) itself. First, we look at the horizontal component 
v* n . It satisfies the following Poisson equation in D R+ 2 An '- 

r\ 271 o in 

-A<" = in D R+2ian , -g- = on z = 0, A, < n = on x = a n . 

Moreover, the C 2 ' Q -regularity of u n implies that the boundary value 4>{z) = v™{R-\-2, z) is bounded 
in C 2 [0, A]. Therefore, as we did with u™, we represent v n n {x,z) = (j)(z)g(x) + v™(x,z) with the 
cut-off function g(x) as in (|3.17|) . The function v™ satisfies 

-A<" = In ■= -Mz)g(x) - <t>{z)g"{x) + in D R+2>an , 

with an exponentially decaying function /„, as follows from (|3.25|) . The boundary conditions are 

dv in 

—p- = on z = 0, A, = on x = R + 2. a n . 
oz 
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The same argument as we used to obtain (|3,21|1 implies that 

ran /"A r 

/ / \v™\ 2 e 2Px dxdz < C / \f n (x,z)\ 2 e 2f3x dxdz < C (3.26) 

JR+2 JO J 

with a sufficiently small < (3 < a. Therefore, in the same vein as we have obtained (|3,24|) and 
()3.25j) . we conclude that 

\v™(x, z)\ < Ce~ ax for x > R + 5 (3.27) 

and 

\Vv™(x, z)\ < Ce~ ax foix>R + 5. (3.28) 

The uniform bound on w 1 ^ now follows, as it satisfies the Dirichlet boundary condition (x , 0) = 

dw m dv m 

w^ix, A) = and the derivative - n = — —2- is exponentially decaying (13.281) . We infer that 

oz ox 

|< n (x, z)\ < Ce~ ax for x > R + 5. (3.29) 

Now we bound u" 11 *. The corresponding vorticity satisfies 

a, .out 

~ Cn ~Bx~ + Un ' VuJ ° nUt ~ aAi0 ° nUt = p & ' V±T ^ in D ^ u = on D *»- 

However, 

|VT° U *| < e (3.30) 
by construction, hence the maximum principle implies that 

\u™\x, z)\ < epq(z) < Ce. (3.31) 

Here the non-negative function q(z) satisfies the boundary value problem 

-aq"{z) + w n {z)q'{z) = 1, q(0) = q(X) = 0. 

We infer from the standard local elliptic estimates up to the boundary, (|3.31|) and (|3.3()|) that 

\Vu; n ut (x,z)\ < Ce in D an (3.32) 

as well. The vertical flow component satisfies 

f),,out f)„..out(_i_ \ 

^ = ^X~ [n D - <U ' (X ' 0) = = °. Q X = °- 

Therefore, the maximum principle implies once again that 

\w° n ut (x,z)\<^fz(\-z)<Ce in D an . 
Hence, the same local elliptic regularity results allow us to conclude that 

\Vw° n u \x,z)\ <Ce in D an . 
In order to bound the horizontal flow component and conclude the proof of Lemma 13.51 we 

f)„.out a„ .out f)„.out 

observe that = - u° ut so that 

Oz oz " oz 

mean-zero condition 



< Ce in D an . However, v™ 1 also satisfies the 



/ v° ut (x, z) = for all -a n < x < a n . 
Jo 



Hence, we have \v° ut (x, z)\ < Ce in D an , and the proof of Lemma 1531 is now complete. □ 
The next lemma implies that the right limit = 0. 
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Lemma 3.6 The right limit 6 + = 0. 



Proof. Let us choose R independent of n so that c an > sup x>R \v n (x, z)\ for all n. Lemma [3.21 
implies that the speeds c n are uniformly bounded below by a positive constant, thus it follows from 
Lemma 13 . 51 that we can find such R > 0. Then the function 4>(x) = Ae~ ax , with a sufficiently small 
a > 0, satisfies 

+ u n • Vcj) > A<j). 

An argument as in the proof of Lemma 12.21 shows that if A is chosen so that Ae~ aR > 1 then 
T n (x,z) < Ae~ ax on the domain x G [R, a n ]. Therefore, the limit T(x,z) obeys the same bound, 
which in turn implies that 0+ = 0. □ 

Finally, we show that under the additional assumption (|1,1U|) the left limit = 1. This is the 
only place in the proof where assumption Ijl.lOj) is used. 

Lemma 3.7 Let us assume that f(T) < (T — 6o)+/\ 2 . Then the left limit 6- = 1. 
Proof. We note that we have for each iGl 

(M(x) - m(x)) 2 



I 



\VT(x,z)\ z dz > 



X 



with M(x) = max 2 T(x, z) and m(x) = min 2 T(x,z). It follows from the maximum principle that 
the function m(x) is non-increasing. Let us assume that 9- < 9q, then monotonicity of m(x) implies 
that m(x) < 8q for all Then we have 

\VT(x,z)\ 2 dxdz> J(M(x)-m(x)) 2 ^> J (T (x , z) - 9 ) 2 + ^ . 

We also observe that 

J \X7T\ 2 dxdz = J \ T- f(T)dxdz. 



so that 



Hence we obtain using (|1.10|) 



(t - e )%^ >((t- 6 -) HT)dxdz > [ (T(x, z) - e ) 2 dxdz 



A 2 -J \ 2 J JX ' - J v K ' ' + A 2 

However, the left side is smaller than the right side unless T = 6q, the case that we have already 
ruled out. □ 

This finishes the proof of Theorem 11.11 



4 Bounds for the initial value problem 

We consider in this section the solutions of the Cauchy problem with general front-like initial data 
and obtain the uniform bounds on the bulk burning rate and other average quantities stated in 
Theorems 11.21 and Theorem 11.31 We prove the first result, and the proof of the second result is 
presented in Section FOl 
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4.1 Bounds in an arbitrary strip 

We prove in this section Theorem ll.2l Let T(t, x, z), u(t, x, z) be the solution of the Cauchy problem 

T t + u • VT = AT + f(T) (4.1) 
uj + u- Vu - aAu + Vp = pTe 
V • u = 0, 

(4.2) 

with initial data Tq(x, z), uo(x). We assume that there exists R > so that To(x, z) = for x > R 
and Tq(x, z) = 1 for x < —R, and that the initial vorticity is bounded in L 2 : 

\u)q{x , z)\ 2 dxdz < +oo. 

The assumptions on the initial temperature To can be relaxed - it simply has to approach one and 
zero at the two ends of the domain sufficiently fast. 

We recall that the bulk burning rate V(t), the Nusselt number N(t) and the average horizontal 
flow U (t) are defined by 

Vit) = \l * f V(s)ds, V{t) = [ f(T) dx ' h 



t 



i r* 

N(t) = - I N(s)ds, N(t) 



U{t) 



10 

if 



|VT| : 



A 

, dxdz 



Ks)||ooGte- 



(4.3) 
(4.4) 
(4.5) 



The laminar front speed Co is defined as the unique c so that equation 

- c $' = $" + /($), $(-oo) = 1, $(+oo) = 
has a solution < < 1 . We recall the statement of Theorem 11.21 

Theorem 4.1 There exists a constant C > so that under the above assumptions on the initial 
data Tq,uq, the following bounds hold 



c - C[p + p 2 } + o(l) < V{t) < c + C[p + p 2 } + o(l) 

2 

+ o(D 



(4.6) 



N(t)< 
U(t) < Cp[l + p] + o(l) 



CO , ^,_ 2 



Cp+\/^ + C 2 p 



as t — > +oo. 



This theorem shows that the coupled problem (|4,1[) is in a sense a regular perturbation of the single 
reaction-diffusion equation with p = 0. The lower bound in (|4,6|) is of interest only for small p when 
the left side is positive. 

Proof. First, we prove the following bounds on N(t) and V(t) in terms of U(t). 



Lemma 4.2 There exists a constant Cq that depends only on the initial data T$ so that 



N{t) < ±V(t) + U(t) + Co 



1 1 

t + Vt 



and 



V(t) < c + U(t) + C 



1 1 



(4.7) 
(4.8) 
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Proof. Define g(T) = T(l — T) and its integral 

dxdz 



R{t) = J g(T)'- 



The idea of using a concave function g(T) in a related context is due to B. Winn 38 . We observe 
that 

f = J m AT d -^ + J V(T)/m^2- / >c0|v3l'^-ni) (4.9) 

with the burning rate V(f) defined in (|4.3|) . Thus 

^ + K(t) > 2 y |VT| 2 — = 2iV(i), 
which after averaging in time becomes 

+ V(t)>2N(t). (4.10) 

In order to obtain an upper bound for the potentially small term R(t)/t in (|4.1U|) we construct 
sub- and super-solutions for T(t, x, z). This construction follows jlJJj. We look for a sub-solution for 
T of the form 

ipi(t,x,z) = ®o(x - c t + x 1 + £i(i)) - qt(t,x,z). 

Here <l>o is the traveling wave in the absence of convection, at p = 0, normalized so that ^(O) = /Jo- 
lt is the unique solution of 

-co$ = *o + /($o), *o(0) = 0o, *o(-oo) = 1, $ (+oo) = 0. 
The functions £i(i) and gi(t, x, z) are to be chosen. In order for ipi to be a sub-solution we need 

G[tk] = ^ + u • - - /(^) < 0. 

We have 

G[M = L% +u& -^-u- Vqi + Agi + /($ ) - /($o - ffi). 

With an appropriate choice of x%, that is, by shifting $o sufficiently to the left we can ensure that 
T (x,z) >&o(x)-q 10 (x) withO<fto(x) < (l-0 o )/2 and q 10 (x) € Then we choose q\{t, x, z) 

to be the solution of 

^ + u • Vgi = Agi, gi(0, x, 2) = qio(x), = at z = 0, A. (4.11) 

The following lemma first proved in ^B] provides a uniform L 1 — L°° decay estimate for q\ that is 
independent of the advection term. 

Lemma 4.3 There exists a constant C > i/iai is independent of the (incompressible) flow u so 
that 

C 

l|9l(*)l|oo < ^lklollLi(D) (4-12) 

/or i > 1. 
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As mentioned above, the main point of the above result is the independence of the constant in (|4.12j) 
from the flow u. We also note that this L 1 — L°° estimate behaves in a one-dimensional way for large 
times, as one would expect for a strip. The factor of A in the denominator is compensated by the 
fact that the L 1 -norm is taken over the strip and not only in x. We postpone the proof of Lemma 
14.31 till the end of this section. 

We can find 5 > so that if $ 6 (1 - 6, 1) and 91 G (0, (1 - o )/2) then /(*<,) < /(*o - <*)• 
Hence we have in this range of <£o : 

G[^\ < + v& . (4.13) 

Furthermore, if <£o G (0) then /(<3?o) = f{&o — 5) = and hence in this range of $0 we have (|4.13|l 
with the equality sign. Finally, if 3> & ($, 1 — &) then \f{$o) - /($o - q)\ < K\q\ and % < -(3. 
Hence G[ipi] < everywhere provided that 



Ut)> |Nt)|U + *''^"°° . (4.14) 
Thus choose 

£ 1 {t) = U(t)t + CVi. (4.15) 
Therefore we obtain a lower bound for T: 

T(t, x, z) > <£ (x - c t + U(t)t + CVi) - qt(t, x, z). (4.16) 

In order to obtain an upper bound we set ip u = &o(x — cot — x<i — ^(t)) + q2(t,x,z) and look for 
£,2(t) and q2(t,x, z) so that G[ip u ] > 0. The constant X2 is chosen so that 

T (x,z) < $ (x - x 2 ) +q2{0,x,z) 

with 92(0, x,z) £ ^(D) and < 92(0, x, z) < 6o/2, as with qi(0,x, z). The function q2(t,x,z) is 
then chosen to satisfy the same advection-diffusion equation (|4.11|) similarly to q\. Hence it obeys 
the same time decay bounds as q%. With the above choice of q2 we have 

G(V>u) = -6*0 + v&o + /( $ o) - /($o + 92). 

Once again, we consider three regions of values for &q. First, if 1 — 5 < $0 < 1 with a sufficiently 
small 6 > then /($o) — f(^o + Q2) > 0, as q2 > 0. Hence G[ip u ] > in this region provided that 
£2 > 0. Second, as 92 < 9q/2 we have /(3>o) = /(^o + Q2) = if < $0 < <5 with a sufficiently small 
5 > 0. Hence G[?/> M ] > in that region under the same condition £2 > 0. Finally, if $0 £ 1 ~~ <5) 
then $' < -/3 with /? > and |/($o) - /($o + 92)1 < ^ll^lloo- That means that G[ip u ] > if we 
choose £2 so that 

^H<?2||oo 



6>Kt)||oo+ , . 

Therefore we choose 

£ 2 (t) = U(t)t + CVt, 
as with £i(t). Therefore we obtain upper and lower bounds 

<S>o(x ~ cot + ft(t) + xi) - ft (t, x, z) < T(t, x, z) < $ (x - c t - - x 2 ) + <ft(t, x, z) (4.17) 

that imply in particular that 

$o(x ~ cot + U(t)t + C [l + Vi])-^< T(t, x, z) < $ (x - c t - U(t)t - C [l + V~t]) + Q (4.18) 
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with a constant Co determined by the initial conditions. Hence, using (|4.17jl - (|4.1H|) and the L 1 - 
bounds ||9j(i)||z,i(D) — Co, j = 1,2, we obtain 



R(t) = / T(l - T) 



dxdz 



CQt-&(t)-X2 



m , ^.dzdx 

r(i-r)— — + 

A 



co*+£i(*)+zi 
co*-^2(*)-a;2 



T(l - T) 



dzdx 



+ / T(1-T)—-<C + (l-*o)dx + (e 1 (t)+6(*))+ / $o(z)^ 

•/cot+£i(t)+xi J A J-oo JO 

< Cb(l + Vt)+2tU(t). 

This together with (|4,1U|) implies that 



+ 2tf(t) + C 



1 1 



> 2N(t) 



(4.19) 



so that (|4~T|) holds. 

Moreover, we have 



P(t) 



1 



f(T(s,x,z))^p- J ds = j 



[T(t,x,z) - T (a:,z)] 



dxdz 



T t (s,x,z] 



dxdz 



ds 



(4.20) 



< i y [$ (z - cot - x 2 - &(t)) + 92(2, a?, z) - ^0(2; + xi) + gi(0, x, z 
Co 1 [ c o t +i2{t)+x 



dxdz 



1 - $ (^))dx + - / §o{x)dx < c + I7(t) + C 



1 1 

Vt + t 



as follows from (|4.17jl - (|4.18jl . This proves Q4.8J) and finishes the proof of Lemma 14.21 □ 
On the other hand we have the following upper bound for U(t) in terms of N(t). 

Lemma 4.4 There exists a constant C > so that for all t > the following inequality holds 



U(t) < C 

Proof. We multiply the vorticity equation 

8uj 



P\/N(t) + — Hwoll^ 



(4.21) 



— + u • Vw - auj = p(e ■ V^T) 



by co and integrate: 

— — J \oj(t, x, z)\ 2 dxdz + a J \ S/uj(t, x, z)\ 2 dxdz = p J u(t, x, z)(e^ ■ X7T)d. 

with e 1 - = (e2, —ex)- 

The Poincare inequality for u(t, x, z) implies then that 



xdz, 



(4.22) 



1 d 
2~dt 



a 



u(t,x,z)\ dxdz + — / \Vuj(t,x,z)\ dxdz<Cp / \VT(t,x,z)\ dxdz. 



Integrating this equation in time we conclude that 



I Voj(s, x, z)\ 2 dxdzds < C 



p 2 7V(t) + -|| Wo |li2 



(4.23) 
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However, as in the proof of Lemma 12.41 we have ||f (i) (i?) — C||^ aJ (^)llL 2 (D)- This, together with 
(1031) implies (|OT|) . □ 

Putting the bounds (|4,21|) and (|4,7j) - (|4,8[) together and using the Cauchy- Schwartz inequality we 
arrive at 



1 1 

t + 7t 



+ 



2N(t) < c + C 
Hence we obtain an upper bound 

N(t)< 

This, together with (|4.21|) implies that 



Cp 



I ^/N{s)ds < c + C 
Jo 



1 1 



+ Cp\N(t). 



+ 0(1). 



(4.24) 



It follows then from (|4.2U|) that 



U(t) < Cp[l + p]+o(l) 



V(t)<c + Cp[l + p]+o(l). 



The lower bound on V{t) in (|4.6|) is proved similarly. This finishes the proof of Theorem 14.11 It 
remains only to prove Lemma l4.HI □ 

Proof of Lemma 14.31 We will show that there exists a universal constant C > so that the 
solution of 



+ u • Vtp = aAifj 



(4.25) 



rp(0,x,z) = ip (x,z) > 0, 
with the Neumann boundary conditions at z = and z = X, and u sufficiently regular, satisfies 

\Mt)h°°(D) <Cn 2 (t)\\^ \\ L1{D) . (4.26) 
Here n{t) is the unique solution of 

n\t) C 



l + n 3 (t)A 3 a\H' 
We multiply (|4.25|) by if) and integrate over the domain D to obtain 

1 d 



2 dt 



(4.27) 



(4.28) 



We now prove the following version of the Nash inequality _29 for a strip of width A in two dimensions: 



\\vnl>c 



(4.29) 



li + A 3 |miiimi 2 

The proof of ()4.29|) is similar to that of the usual Nash inequality. We represent ip in terms of its 
Fourier series-integral: 

dk 



n>Cr M 7F 
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where 

4>n(k) = — J e~ lkx cos (^-J ip(x,z)dxdz. 



X 

Therefore we have 



\Mk)\<j\m L i. (4.30) 



The Plancherel formula becomes 

dkdpdxdz 



[ mx,z)\ 2 dxdz= f e ik ^co S (^)cos(^)Uk)^ m (py- 



n,m>0 
n>0 



(2t) 5 



and similarly 

n>0 " 



Let p > be a positive number to be chosen later. Then using the above Plancherel formula we 
write 

U\\l=i + n, 

with the first term that is bounded using (|4.«S0j) 

A [ „.m2<» „ C\p([\ P \ + i) Op(X P + i) 



, |t|</i ^VT A - A 

0<n<pA" I ^ I — ^ 

The rest is bounded by 

CX^ f ( ' , 47r 2 n 2 \ , 7 ,,,,, „ C 



P tXnJken V a / P 



n>0 

Therefore we have for all p > 0: 

^ Cp(Ap+l) j+ C , 
A ,cr 

We choose p so that 

AIIV^IH 



P i i i 2 



II 

and obtain 



2 < C||V^||f / A 4 / 3 ||VV>||f \ 2 C7||V^|| 



lll^llf 



M 3 \M\l /S V IMlf / A2/3|| V ^||f 

^_ii.,.iiJ/3 ||v^||f + CA 2 / 3 ||W|| 4/3|I - I|2/3 



A 2 /3 



1 ' 



This is a quadratic inequality ax 2 + bx — c > with x = WVipW^ 3 , a = CX 2 ^ 3 \\ip\\^ 3 , b = — ^ • ' 
and c = HV'lli an d hence 

-6 + Vb 2 + 4ac 2c c 

# > = ; > 



2a b + V6 2 + 4ac Vfr 2 + 4ac' 
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This implies that 



and therefore 



8/3 \ — 3/2 ^ —1 

II Will > cii^m f + 4A 2 / 3 ||^||f Ml) > c|« (Mi + aimumiI 



> CA 2 ||^II| 



Hence (j4"2U)) indeed holds. 

We insert (|4.29|) into (|4.28j) and use the conservation of the L 1 -norm of tp (recall that the initial 
data is non-negative) obtain 

> < - 4 C ^"| (4.3!) 



ll^olll + A 3 ||Vnl ! 1 1 '■ ' 1 1 2 

Integrating (|4.31f) in time we have 



4 A 3 ||i/>olli , 1 



zHt) 



It Mia ~z{t) 

where z(t) = \\ip(t)\\2/\\ipo\\i, anci thus 

gHO < 1 (432) 
1 + A 3 z 3 (i) - CaAV 1 j 

The function on the left side of ()4.32j) is monotonically increasing and hence we have 

| 2 <n(i)||Vo||i, (4.33) 



where n{t) is the solution of (|4.27|) . 

Let us denote by Vt the solution operator for (|4.25|) : = Vt^Po- Then ()4.33|) implies that 

||'Pi||xi_ + x2 < n(t). The adjoint operator VI is the solution operator for 

^ - u • V-0 = ctAt/- (4.34) 
V>(0, x) = tpo(x), x <E D 

with the Neumann boundary conditions at z = 0, A. Note that the preceding estimates rely only on 
the anti-symmetry of the convection operator u- V. Therefore we have the bound \\Vf Wli^l? < n(t) 
and hence ||7^t ||z, 2 ^z, 00 ^ n (t) so that 



< n(t/2)U(t/2)\\ L 2 <n 2 (t/2)\\^\ 



and thus (|4.26j) indeed holds. 

The estimate (|4.12j) follows from the observation that for large i > 1, when n(t) is small, we 
have the bound 

n 2 (t)< ( ' 



This finishes the proof of Lemma 14.31 □ 
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4.2 Bounds on the burning rate in a narrow domain 

We recall that no non-planar traveling fronts do exist in the reactive Boussinesq problem in a narrow 
vertical strip when gravity is sufficiently small |32l 1331 ITl] . Moreover, solutions with general front- 
like initial data become asymptotically planar in the long time limit |14j . We extend now this result 
to the inclined cylinders. More precisely, we have the following result (this is a re-statement of 
Theorem [Oil. 

Theorem 4.5 Let e = (ei, e%) be the unit vector in the direction of gravity and let pj = pej, j = 1, 2. 
There exist two constants Xq and pq so that if the domain is sufficiently narrow: A < Ao and gravity 
is sufficiently small: p < po then the burning rate is bounded by 

V(t) < c + Cp 2 + o(l) as t -> +oo. (4.35) 

Moreover, the front is nearly planar in the sense that 



1 



N z (t) = - / \\T z (s)\\ids < Cpi + o(l) as t -» +oo. (4.36) 
t Jo 

The key point in Theorem l4.5l is that the bounds in (|4.35|) and (|4.36jl are independent of the gravity 
strength p\ in the direction parallel to the cylinder. 

Proof. Multiplying the vorticity equation by oj and integrating by parts we obtain 



\d_ 
2~dt 



J \u\ 2 dxdz + a j \V uj\ 2 dxdz = p2 j T x uodxdz — p\ J T z udxdz. (4.37) 



The Poincare inequality applies to uj(x,z) with the Poincare constant proportional to 1/A. Hence, 
if A < Aq and pj < po, (|4.37l) implies that 



1 d 
2~dt 



j \uj\ 2 dxdz + ^ j \Vu\ 2 dxdz <Cp\ J \T X \ 2 dxdz + C p\ j \T z \ 2 dxdz. (4.38) 



We now differentiate the equation for T(t, x, z) in z to get 

BT 

+ u • VT Z + u z ■ VT = AT Z + f'(T)T z . 

at 

Multiplying this equation by T z we obtain 

\jt j \ Tz \ 2dxdz + f NT z \ 2 dxdz + j T z u z -VTdxdz = j f(T)T 2 < M j T 2 dxdz. (4.39) 
The last integral on the left side is bounded by 

Tu z ■ X7T z dxdz 



T z u z • VTdxdz 



< 2 / \u z \ 2 dxdz + — \VT z \ 2 dxdz. 



This, together with incompressibility of u, the Poincare inequality for T z and Q4.39JI imply that 

\~d\J \ T z\ 2(ixdz + \ J \ VT z \ 2 dxdz < 4 j \uj\ 2 dxdz, (4.40) 



provided that A < Ao- Combining (|4.38|) and Q4.4UJI and using the Poincare inequality for u and T z 

P II 2 - 

[ ^ll2- 

idn c 



once again, we obtain the following inequalities for U(t) = \\u!(t)\\?, and N z (t) = ||T 2 " 2 - 



., (U +^<Cp 2 N x (t) + Cp 2 N z (t) 



31 



and 



ldN z (7 

+ TnN z < 40. 



2 dt A 2 
Hence, the function Q = N z + satisfies 



2~dt 



+ 



C_ 

A 2 



4-Cp? 



Q < Cp 2 2 N x (t). 



Therefore, we have 



Q(t) < Qoe" 7 * + Cp\ [ e^^N x (s)ds 
Jo 

with 7 > provided that C / A 2 > 5 and Cp\ < 1 . We conclude that 



Q(t) :-- 



t 



Q{r)dT < 



Qc 



-'JT 



dr + 



( FT 



o J o 



< + £d t 
- t t 7 



^ + C> 2 
* 7 



A", (.s)c / (-~ r <h,ls < + 1-L±x t (t) < Cp\ + -j-. 



The last inequality above follows from the bound on N(t) in Theorem 14.11 Now, the bound (|4.36|l 
in Theorem 14.51 follows. Then, H4.38j) together with H4.36() and the same uniform bound on N(t) in 
Theorem 14. II imply that 

j£\\Vu(a)\\ld8<Cf% + Q 
We recall that \\v (i)llL°°(D) ^ C|l^ 7a; llL 2 (D) - this, together with the above, imply that 



U(t)<Cp 2 + ^. 



(4.41) 



Finally, using (|4~2U)) and (jOT|) we obtain (|Q5|) . □ 



References 

[1] M. Abel, A. Celani, D. Vergni and A. Vulpiani, Front propagation in laminar flows, Physical 
Review E, 64 6307 (2001). 

[2] M. Abel, M. Cencini, D. Vergni and A. Vulpiani, Front speed enhancement in cellular flows, 
Chaos 12, p. 481. 

[3] B. Audoly, H. Berestycki and Y. Pomeau, Reaction diffusion en ecoulement stationnaire rapide, 
C.R.Acad. Sci., Ser. IIB, 328, 255-262. 

[4] M. Belk, B. Kazmierczak, V. Volpert Existence of reaction-diffusion-convection waves in un- 
bounded cylinders, Preprint, 2004. 

[5] H. Berestycki, The influence of advection on the propagation of fronts in reaction-diffusion 
equations, in Nonlinear PDEs in Condensed Matter and Reactive Flows, NATO Science Series 
C, 569, H. Berestycki and Y. Pomeau eds, Kluwer, Doordrecht, 2003. 

[6] H. Berestycki and F. Hamel, Front propagation in periodic excitable media, Comm. Pure Appl. 
Math. 55, 2002, 949-1032. 



32 



[7] H. Berestycki, F. Hamel and N. Nadir ashvili, The speed of propagation for KPP type problems 
in periodic and general domains, Preprint, 2003. 

[8] H. Berestycki, F. Hamel and N. Nadirashvili, Elliptic eigenvalue problems with large drift and 
applications to nonlinear propagation phenomena, Preprint, 2003. 

[9] H. Berestycki, B. Larrouturou and L. Nirenberg, A nonlinear elliptic problem describing the 
propagation of a curved premixed flame, in Mathematical Modeling in Combustion and Related 
Topics, C.-M. Brauner and C. Schmidt-Laine, eds., NATO ASI Series, Kluwer, 1988. 

[10] H. Berestycki, B. Larrouturou and P. L. Lions, Multi-dimensional traveling wave solutions of a 
flame propagation model, Arch. Rational Mech. Anal., Ill, 1990, 33-49. 

[11] H. Berestycki, B. Nicolaenco and B. Scheurer, Traveling wave solutions to combustion models 
and their singular limits, SIAM Jour. math. Anal., 16, 1983, 1207-1242. 

[12] H. Berestycki and L. Nirenberg, Traveling fronts in cylinders, Annales de 1'IHP, Analyse non 
lineare, 9, 1992, 497-572. 

[13] P. Constantin, A. Kiselev, A. Oberman, L. Ryzhik, Bulk burning rate in passive-reactive diffu- 
sion, Arch. Rat. Mech. Anal. 154, 2000, 53-91. 

[14] P. Constantin, A. Kiselev and L. Ryzhik, Fronts in reactive convection: bounds, stability and 
instability, Comm. Pure Appl. Math., 56, 2003, 1781-1803. 

[15] A. de Wit, Fingering of chemical fronts in porous media, Phys. Rev. Let., 87, 054502. 

[16] A. Fannjiang, A. Kiselev and L. Ryzhik, Unpublished notes, 2002. 

[17] R. Fisher, The wave of advance of advantageous genes, Ann. Eugenics, 7, 1937, 355-369. 

[18] M. Freidlin and J. Gartner, On the propagation of concentration waves in periodic and random 
media, Soviet Math. Dokl., 20, 1979, 1282-1286. 

[19] M. Freidlin, Geometric optics approach to reaction-diffusion equations, SIAM J. Appl. Math., 
46, 1986, 222-232. 

[20] D. Gilbarg and. Trudinger, Elliptic Partial Differential Equations of Second Order, Springer- 
Verlag, 1983. 

[21] F. Hamel, Formules min-max pour les vitesses d'ondes progressives multidimensionnelles, Ann. 
Fac. Sci. Toulouse Math., Serie 6, 8, 1999, 259-280. 

[22] S. Heinze, G. Papanicolau and A. Stevens, Variational principles for propagation speeds in 
inhomogeneous media, SIAM J. Appl. Math. 62, 2001, 129-148. 

[23] L. Kagan and G. Sivashinsky, Flame propagation and extinction in large-scale vortical flows, 
Combust. Flame 120, 2000, 222-232. 

[24] L. Kagan, P.D. Ronney and G. Sivashinsky, Activation energy effect on flame propagation in 
large-scale vortical flows, Combust. Theory Modelling 6, 2002, 479-485. 

[25] A. Kiselev and L. Ryzhik, Enhancement of the travelling front speeds in reaction-diffusion 
equations with advection, Ann. Inst. H. Poincare Anal. Non Lineaire 18, 2001, 309-358. 



33 



[26] A.N. Kolmogorov, I.G. Petrovskii and N.S. Piskunov, Etude de l'equation de la chaleurde 
matiere et son application a un probleme biologique, Bull. Moskov. Gos. Univ. Mat. Mekh. 1 
(1937), 1-25. (see [30] pp. 105-130 for an English transl.) 

[27] A. Majda and P. Souganidis, Large scale front dynamics for turbulent reaction-diffusion equa- 
tions with separated velocity scales, Nonlinearity, 7, 1994, 1-30. 

[28] S. Malham and J. Xin, Global solutions to a reactive Boussinesq system with front data on an 
infinite domain, Comm. Math. Phys., 193, 1998, 287-316. 

[29] J. Nash, Continuity of solutions of parabolic and elliptic equations, Amer. Jour. Math., 80, 
1958, 931-954. 

[30] Dynamics of curved fronts, P. Pelce, Ed., Academic Press, 1988. 

[31] G. Papanicolaou and X. Xin, Reaction diffusion fronts in periodically layered media, Jour. Stat. 
Phys., 63, 1991, 915-932. 

[32] R. Texier-Picard and V. Volpert, Problemes de reaction-diffusion-convection dans des cylindres 
non bornes, C. R. Acad. Sci. Paris Sr. I Math. 333, 2001, 1077-1082 

[33] R. Texier-Picard and V. Volpert, Reaction-diffusion-convection problems in unbounded cylin- 
ders, Revista Matematica Complutense, 16, 2003, ... 

[34] N. Vladimirova, R. Rosner, Model flames in the Boussinesq limit: the effects of feedback, Phys. 
Rev. E., 67, 2003, 066305. 

[35] V.A. Volpert and A.I. Volpert, Location of spectrum and stability of solutions for monotone 
parabolic system, Adv. Diff. Eq., 2, 1997, 811-830. 

[36] V.A. Volpert and A.I. Volpert, Existence and stability of multidimensional travelling waves in 
the monostable case, Israel Jour. Math., 110, 1999, 269-292. 

[37] V.A. Volpert and A.I. Volpert, Spectrum of elliptic operators and stability of travelling waves, 
Asymptotic Anal., 23, 2000, 111-134. 

[38] B. Win, Ph.D. thesis, University of Chicago, 2004. 

[39] J. Xin, Existence of planar flame fronts in convective-diffusive periodic media, Arch. Rat. Mech. 
Anal., 121, 1992, 205-233. 

[40] J. Xin, Existence and nonexistence of travelling waves and reaction-diffusion front propagation 
in periodic media, Jour. Stat. Phys., 73, 1993, 893-926. 

[41] J. Xin, Analysis and modelling of front propagation in heterogeneous media, SIAM Rev., 42, 
2000, 161-230. 

[42] Ya.B. Zeldovich, G.I. Barenblatt, V.B. Librovich and G.M. Makhviladze, The Mathematical 
Theory of Combustion and Explosions, Translated from the Russian by Donald H. McNeill. 
Consultants Bureau [Plenum], New York, 1985. 



34 



